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ABSTRACT 


Flow-induced vibration (FIV) on two tandem cylinders with forced convection is numerically investigated at a constant Re = 150. Elastically 
mounted cylinder with four different values of corner radii (r* = r/R: r = radius of fillet; R = radius of circle) = 0 (square cylinder), 0.25, 
0.75, and 1 (circular cylinder) with two spacing ratio (L/D) = 4 and 2 is studied. Transverse oscillations are generated from the cylinder 
having non-dimensional mass (m*) = 10. The structural damping coefficient is assigned a zero value with varying reduced velocity 
Ur = 2 — 10. The two-dimensional incompressible Navier-Stokes and energy equations are solved together with Newton’s second law gov- 
erning the motion of the cylinders. Both cylinders’ surfaces are maintained at a higher constant temperature of T* = 1, and incoming flow is 
set to be at T* = 0 with Prandtl number (Pr) = 0.7. The effect of r* and L/D is observed on the flow structure and FIV parameters. Flow 
characteristics at L/D = 4 such as steady flow, reattachment, and unsteady flow are examined. A “shift” in vibrational amplitude is noted 
from r* = 1 and 0.75 to r* = 0 and 0.5, respectively. The downstream cylinder (Dc) experiences a hike in vibration amplitude due to the 
impingement of vortex shedding from the upstream cylinder (Uc). r* = has 18.1% higher vibrational amplitude than r* = 0 at their respec- 
tive lock-in regimes for Dc. For L/D = 2, vortices from upstream and downstream cylinders interact to form C(2S) and 2S types of vortex 
shedding. Different regimes, such as single body, reattachment, and co-shedding, have been observed while changing L/D. r = 0.75 results 
in 13.3% higher oscillation amplitude as compared to r. = 0.5 for Dc. The average Nusselt number (Nag) strongly depends on flow topol- 
ogy, corner radius, and vibrational amplitude A/D. At low L/D, heat transfer from the downstream cylinder is plummeted due to rolling of 
shear layers over the cylinder. There is a significant change in Nuag due to higher vibration; for example, increase in 10.71% change is 
observed from Ur = 2 to Ur = 6 for r = and L/D = 4. Corner radii also alter the Nulayg as a decrease in 27.39% from r* = 1 to r* = 0 at 
Ur = 10 and L/D = 4 (Uc). 
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NOMENCLATURE F, Normalized natural frequency 
k Spring constant (N/m) 
A Maximum amplitude (m) m Mass (kg) 
c Damping coefficient (N s/m) n Mass ratio 
Cp specific heat capacity (J/kg K) Mic Non-heated cylinder 
Ct Coefficient of lift Nu Nusselt number 
Cp Average coefficient of drag Nuag Average Nusselt number 
C™ Maximum coefficient of drag Pressure (Pa) 
D Diameter of cylinder (m) Pr Prandtl number 
De Downstream cylinder (m) r= Corner radii (r is the radius of the corner, R is the radius 
nod Modulation frequency (s J of the Cylinder) 
fa Natural frequency (s ) Re Reynolds number 
5. Frequency of oscillation E5 Ri Richardson number 
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Ts Surface temperature (K) 
Tx, Flow temperature (K) 
St Strouhal number 
„x. velocity 
Uc Upstream cylinder 
Ux Free stream velocity (m/s) 
Ur educed velocity 
v y-velocity 
F. direction amplitude 
Thermal conductivity (W/m K) 


1 Viscosity (Pa s) 

€ Damping ratio 

p Density (kg/m*) 
l. INTRODUCTION 


Flow-induced vibration (FIV) is found to affect the structural 
integrity of marine structures (Xu et al, 2022), heat exchangers 
(Cheng et al., 2009; Ji et al, 2022; and Qiu et al., 2021), ocean pipelines 
(Ali et al, 2022b; Gray et al, 2015; and Rao et al, 2017), and transmis- 
sion lines (Gao et al., 2021; Zhu et al., 2019). FIV on an isolated bluff 
body, for example, single cylinder, has been investigated in existing 
kinds of literature (Williamson and Govardhan, 2004; Jauvtis and 
Williamson, 2003). For a single elastically mounted cylinder, its 
motion under a cross-flow configuration is governed by various factors 
including spring constant (k), damping constant (c), and natural fre- 
quency (fa) (Chaplin et al., 2005; Liu et al., 2020). For a cylinder in 
cross-flow configuration, vortex shedding can generate vibrations if Re 
is not very low (Wang et al., 2017; Xiao et al., 2011). Williamson and 
Roshko (1988) divided vortex shedding type on the basis of 
amplitude-phase graph. High-amplitude vibration occurs when vortex 
shedding frequency (f;) approaches fn. This phenomenon is known as 
“lock-in synchronization” or “frequency synchronization” (Nepali 
et al., 2020; Sarpkaya, 2004). Khalak and Williamson (1997) studied 
the effect of combined mass-damping (m*€) on vortex shedding and 
amplitude modes. In the amplitude-phase graph, higher mass ratio is 
associated with the “initial” (2S vortex shedding mode, two singlets in 
a single period) and “lower” (2P vortex shedding mode, two pairs in 
single period) branches. However, in low mass ratio, lock-in region 
increases, and an “upper” branch appears. Low mass ratio leads to a 
broader lock-in region with only one “upper branch” of 2P vortex 
shedding mode. 

FIV on two bluff bodies, for example, tandem cylinders, is a com- 
plex phenomenon. It is of course physically more interesting than FIV 
on a single cylinder. For tandem cylinders, the motions of both cylin- 
ders are intimately coupled through the hydrodynamic forces 
(Borazjani and Sotiropoulos, 2009; Mussa et al., 2009). FIV on two 
tandem cylinders can be classified into different regimes based on the 
vibration amplitude and flow pattern observed. Systematic study of 
flow around tandem cylinders gives deep insight into flow structure. 
Different in spacing ratio L/D (L: center to center distance; D: diame- 
ter) results in different vibration regimes: single-body regime, reattach- 
ment regime, and co-shedding regime (Zdravkovich, 1996; 1987; Zhou 
and Mahbub Alam, 2016). For/D < 2.5, single-body regime exists 
when both cylinders act as if they are an isolated bluff body in flow. 
For 2.5 < L/D < 5, flow separation occurring from upstream cylin- 
der (Uc) reattaches on the downstream cylinder (Dc), referred to as 
reattachment regime. For large spacing L/D > 5, vortices are shed 
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independently from each cylinder and are referred to as the co- 
shedding regime (Khan ef al., 2022; Sumner, 2010). For large spacing 
ratio, wake galloping and hitting of vortices generated from the Uc on 
the Dc are observed (Griffith et al., 2017; Prasanth and Mittal, 2009). 
Lin et al (2014) numerically investigated the effect of varying 
m* = 0.5 — 4 and Re=5-100 on isolated and tandem cylinders. In 
addition, the symmetrical vibration regime, an exceptional regime, is 
observed for tandem cylinders. Asymmetrical regime has been discov- 
ered even at low Re = 10 and Ur > 7. In the asymmetrical regime, the 
oscillation axes shift from the initial position and a “shift” occurs. 

The above-cited studies were performed on circular cylinders, 
which is a special case of square cylinder with filleted corner radius 
exactly equal to half of its side length. More generally, for square cylin- 
ders, the sharp corners can be filleted to a specific corner radius to con- 
trol flow separation location with the aim to alter the hydrodynamic 
forces and therefore mitigate FIV. For a square cylinder, flow separa- 
tion generally occurs at the leading and trailing edges (Hu et al., 2006; 
Kumar and Sen, 2021; and Sarout et al., 2022b). As the sharp corners 
are rounded off, flow separation location changes due to the Coanda 
effect. For flow around a single square cylinder, Ajith Kumar ef al. 
(2009) conducted a study on the influence of corner radius experimen- 
tally. The experiment was performed on six different r = 0.0-0.5 
(where r* = r/R with corner radius r and half side length of a square 
cylinder R) in a water channel at Re=5200. Sarout ef al. (2022a) 
numerically investigated the effect of ? = 0-1 and Re = 100, 150, and 
200 on flow characteristics. Four different modes were observed; low- 
amplitude mode, beating mode, lock-in mode, and monotonic oscilla- 
tion mode. The effect of r* manifested clearly in the shedding regime 
and shedding cycle timing. Zhao and Zhao (2019) numerically exam- 
ined the effect of r* within the range of 0 to 0.5 on vibration modes at 
a constant Re=200 where galloping and lock-in were observed. 
Galloping was completely suppressed for r = 0.1. Amplitude for lock- 
in increased due to a change in r* from 0 and 0.5. In further analysis, 
beating was divided into two parts: type A and type B. In another study, 
Miran and Sohn (2015) numerically investigated the effect of corner 
radius of a square cylinder at Re = 500 for transient two-dimensional 
(2D) laminar with turbulence treated via large-eddy simulations. Six 
different r = 0.0-0.5 were examined for its effect on St. For r within 
the range of 0 to 0.2, a huge leap in St is registered and then gradually 
rose for the range of 0.2 < r* < 0.5. Sen and Mittal (2011) numerically 
investigated flow around a square cylinder with for Re in the range of 
60 < Re < 250. Three main regimes were observed with increase in Re: 
primary lock-in, desynchronization, and secondary lock-in. 2S type of 
vortex shedding was observed in initial branch. In lower branch, both 
2S and C(2S) are detected. Vortex shedding in secondary lower branch 
was through 2S and 2P modes. Zhao et al. (2014) experimentally stud- 
ied the effect of square cylinder angle of attack « = 0°, 20°, and 45° 
with varying U. During the study, the effect of x and Ur on different 
five regimes was observed. For flow around filleted tandem cylinders, 
Adeeb et al. (2018) employed immersed-boundary lattice Boltzmann 
method to address the effect of F in the range of 0.1-0.5 and L/D = 
1.5 — 10 at Re = 100. L/D and r affected both the separation and reat- 
tachment points with changes in flow regimes. 

The interaction of flow with other factors such as heat transfer 
and structural elasticity might be crucial in understanding the FIV 
mechanism and mitigating the damage it causes. The most prevalent 
natural phenomena are heat transfer. A detailed review can be found 
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in Ali et al. (2021); Derakhshandeh and Alam (2019). Flow around 
heated cylinder with FIV demonstrates new findings, which govern 
the nature of temperature variation around the cylinder by Sun et al. 
(2021). Geometrical shape is a major factor in changing topology of 
low around the cylinder and heat transfer phenomenon. Alam et al. 
(2021) numerically investigated the flow heat transfer over a fillet cyl- 
inder at Re = 40-180. For r* = 0.0-1.0, flow topology, aerodynamics 
forces, and heat transfer were explored with four main flow patterns: 
steady trailing-edge separated flow, unsteady trailing-edge separated 
ow, unsteady separation-bubble flow, and unsteady leading-edge sep- 
arated flow. Alam et al. (2021) reported the effect of corner radius and 
the angle of attack & on heat transfer at a constant Re = 150. For 
Corner radius in the range of r* = 0-1, « is varied from 0° to 45°. 
With varying r and 4, four flow patterns were explored: A, B, C, and 
D. Change in heat transfer due to mainly ascribed variations in the 
boundary layer, side-surface bubble, primary wake bubble, and sec- 
ondary wake bubble. Mahir and Altaç (2008) numerically investigated 
the heat transfer in laminar unsteady flow around tandem cylinders at 
e= 100 and 200. For varying L/D = 2 — 10, Nuavg was obtained for 
Uc and De. At higher spacing ratios L/D > 4, Uc acts as an isolated 
cylinder. Dwivedi and Dhiman (2019) numerically inspected heat 
transfer and flow interaction between tandem cylinders at Re= 100 
and Pr = 0.7. Energy equation was solved for five different L/D = 2.5 
—5.5 by using finite volume method. Mean drag coefficient and aver- 
age Nusselt number of Uc was usually higher than De, while for 
higher spacing ratios, change became insignificant. Ali et al. (2022a) 
investigated the forced convection on tandem cylinders at Re = 100 
and varying Ur = 2 — 14 with two degree of freedom oscillating cylin- 
ders. Three different regimes were observed for varying spacing ratio 
L/D = 1.5 — 6: single body, reattachment, and co-shedding regime. 

In the last few years, research on the effect of the corner radius of 
a square cylinder’s FIV has been increased. Existing studies addressing 
the effect of corner radius of square cylinder’s FIV is very limited. 
Most studies are confined to fixed cylinders with constant spacing 
ratio. In addition, coupling to forced convection was not investigated. 
The studies on different cross section, spacing ratio, mounting, and 
heat transfer are summarized in Table 1, with key contribution of past 
researchers. Some of studies have adopted filleted radius along with 
forced convection although the effect of WIV (Wake Induced 
Vibration) and VIV has been omitted. To understand the flow charac- 
teristics and heat transfer, behaviors cannot be predicted on static cyl- 
inders. These features of configurations have not been explored in 
past. A systematic investigation is required on the effect of F and L/D 
on flow characteristics, aerodynamic forces, and heat transfer from 
oscillating filleted cylinders. In present work, the effect of corner radius 
has been studied for tandem square cylinders with a constant 
Re = 150 and varying Ur. Two different values of spacing ratio L/D = 
2 and 4 have been analyzed to understand the physics of tandem cylin- 
ders from r* = 0 (square cylinder) to 1 (circular cylinder). Both cylin- 
ders are maintained at constant wall temperature to examine heat 
transfer along with FIV. 


II. MATHEMATICAL MODELLING 
A. Model configuration and boundary assignments 


Figure 1 shows the schematic diagram of computational domain 
and model configuration. Cylinders are positioned in tandem arrange- 
ment with center-to-center distance (L = 2 and 4), length 2R = D and 
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height 2R = D vertically centered in a large computational domain of 
width S = 70D and height H = 40D constituting a blockage ratio of 
2.5%. Corners have rounded to generate the filleted geometry as 
0< r* < for upstream and downstream cylinder. Both cylinders are 
allowed to vibrate in y direction independently due to hydrodynamic 
forces and constrained in x direction. Mass (m) is actual mass of cylin- 
der, spring constant (k), and damping constant (c) govern the motion 
of cylinders. For heated cases, both cylinders are maintained at con- 
stant temperature T* = 1 and free stream temperature at T* = 0. 
Fluid enters the computational domain at constant density p, viscosity 
v, free stream velocity U, and T,,. The Reynolds number of free 
stream is fixed at Re = 150. 


B. Governing equations 

The incompressible flow and heat transfer around a square cylin- 
der is governed by conservation of mass, momentum, and energy, 
respectively, expressed mathematically as 


=a (1) 
Ox" Oy 
ðu ðu ðu Ldap Pu Pu 
5 ax TY’ Sy 5 dx 05 90 (2a) 
ov av ðv 10p Pv Pv 
att “axe t ay pay 05 ta 5) 
ar, ar, „% (PT PT 60 
a “ax ay Noe o) 


where ii, p, and T are correspondingly velocity, pressure, and tempera- 
ture. For the temperature range considered, fluid properties: density p, 
viscosity 4, specific heat capacity cp, and thermal diffusivity a can be 
assumed constant. Buoyancy effects are not considered as this research 
work is proposed at Ri = 0, which is the ratio between buoyancy force 
and viscous forces. Gourari ef al. (2019) stated buoyancy forces can be 
neglected at Ri < 0.1; hence, forced convection dominates in this 
region. Hence, gravity effect can also be neglected on the basis of Ri 
(Zafar and Alam, 2020). These governing equations are subjected to 
the initial conditions of = 0, v= 0, and T = Tx. The associated 
boundary conditions to the domain are the Dirichlet at the inlet 
(u = Ux, v = 0, T = Te) and Neumann at the outlet (44 = 97 = 0), 
no-slip and no penetration for upper and lower boundaries (u = Ux, 
v = 0, T = T), and the zero output pressure gradient (p = Pout). At 
the cylinder surfaces, no-slip and constant temperature T* = 1 is 
enforced. At low Reynolds number of Re = 150 (Re = 0, the 
flow is reasonably two-dimensional. For the isothermal case and since 
the velocity of the flow is small below the Mach limit of 0.3, the fluid 
indeed is incompressible and the energy equation is uncoupled and 
not involved. However, for the non-isothermal case, the energy 
became coupled, and hence, the density of the flow is followed with 
ideal gas and changes. In general term, the flow appears to be incom- 
pressible but with slight change in density specially near the cylinder 
surface this is not true. Generally, the Mach number is still very small 
to cause any shocking leading to shock wave formation and associated 
density or pressure jump. The flow is considered to be incompressible 
as density and viscosity are constant with temperature, and at 
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TABLE l. Summary of previous studies over different types of cross section with heat transfer. 


Spacing Cross Heat 
Researchers Method ration section Re Mounting transfer Contribution 
Kumar and Sen (2021) Num 5 Square 100 1DOF-1 DOF No Effect of square cross section on 
VIV (Vortex Induced Vivration) 
Khan et al. (2022) Num 1.5, 3.5 Circular 100 1 DOF-1 DOF No Diameter ration and spacing 
and 5.5 ration dependence 
Sun et al. (2021) Num 4 Circular 150 1 DOF-1 DOF Yes Forced convection 
along with FIV 
Zhang et al. (2019) Num 1-8 Filleted 100 Fixed Yes Forced convection 
with rounded corners 
Mittal and Kumar Num 5.5 Circular 100 1 DOF-1 DOF No Beating phenomenon evolution 
(2001) with spacing ration 
Zhu et al. (2019) Num 2, 4, and 6 Circular / Prism 100 Fixed- 2 DOF No Different cross section at upstream 
changes shedding pattern 
Mahir and Alta¢ Num 2, 3, 4, 5, 7, Circular 100 Fixed Yes Effect of spacing on heat transfer 
(2008) and 10. 200 from cylinders 
Jepali et al. (2020) Num 5 Square 40-200 2DOF-2DOF NO VIV of square with an effect of Re 
Adeeb et al. (2018) Num 1.5-10 Filleted 100 Fixed No Effect of corner and spacing 
ration combined 
orazjani Num 17 2, 3 Circular 200 2 DOF DOF No Proximity range explored 
and Sotiropoulos (2009) 
Lin et al. (2014) Num 1.1-4 Circular 5-100 1 DOF-1 DOF No Asymmetric and symmetric 
regimes in FIV 
Alam et al. (2003) Exp 1-8 Circular 65000 Fixed No Fluctuating pressure distribution 
with spacing ratio 
Alam et al. (2016) Exp 4.6-6 Square 13000 Fixed No Staggered angle configuration 
Bao et al. (2012) Num 5 Circular 150 2 DOF-2 DOF No Varying natural frequency 
Chen et al. (2018) Num 1.2-5 Circular 100 1 DOF-1 DOF No Vibration response and 
galloping mechanism 
lu and Zhou (2008) Exp 1.2-2 Circular 7000 Fixed Yes Staggered position with 
heat transfer 
Huera-Huarte and Exp 2-4 Circular ~12000 1DOF-1 DOF No Frequency response 
earman (2011) of experimental vibrations 
Mohany et al. (2014) Exp 1.2-4.5 Circular 10000-300000 Fixed No Flow acoustic in turbulence regime 
himan et al. (2006) Num — Square 1-45 Fixed Yes Effect of varying Re and Pr 
zadpanah et al. (2018) Num = Circular 150 1-DOF Yes Varying Damping ratio along 
with heat transfer 
Sarout et al. (2022a) Num — Filleted 100-200 1 DOF Yes Corner radius with varying 
Re leads to higher heat transfer 
Jwivedi and Dhiman (2019) Num 2.5-5.5 Circular 100 Fixed Yes Thermal buoyancy 
with spacing ration 
e = 150, free stream velocity is relatively low to alter the flow proper- The initial conditions imposed on the square cylinder are Y = 0 and 
ties (Dhiman et al., 2006; Mittal et al., 1997). x =0. 

The upstream and downstream cylinders vibrate transversely The governing equations can be non-dimensionalized by intro- 
under the action of various forces. The motion of the cylinders is gov- ducing dimensionless parameters: f = tU. 72 N, Y*=Y 70 2R), U. 
erned by Newton’s 2nd law of motion as / U, Uu. =, p* p D, and T* = (T — Tu) /(Ts — Too). 

PY dy The dimensionless governing equations are 
mag tea E (4) 
where Y is the y direction displacement of the square cylinder under the aut | Ov" =0, (5) 
influence of the y direction hydrodynamic force per unit length (Fy). Or Oy 
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FIG. 1. (a) Computational domain containing two elastically supported tandem cylinders, (b) pictorial depiction of cylinder changing corner radius from square cylinder to circular 
cylinder with intermediate stages, and (c) meshing of computational domain with zoomed-in. 
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ot* Ox* oy* 5 Re (Ox: Oy* pU2,(2R) 
97. iw aT’ F ƏT* 1 PT PT o These are integrated to obtain the average drag and lift coefficients as 
ðt ð / Pr. Re \ x? * dy? }’ Col 
D (11) 
de Y dy* 2 201 N f 
+ 4nFné (An.)? ; 8 2 N ali 
Be te ge OM e 6) T = Laat, a2) 
where the dimensionless parameters are the reduced velocity Ur = 5 
the Reynolds number Re = 5 the Prandtl number Pr = a damp- C. Solution procedure 
ing ratio & =, mass ratio m* . p* g; Am, normal- The governing equations for the discretized flow domain sub- 
ized saatacal pa uency F, = 7. 5 841 2 ſecuen⸗ jected to appropriate flow conditions were solved using Ansys Fluent, 
1 ene Ee P ERGY which is a finite volume-based solver. The transient, convective, and 
In 2. mw diffusive terms are discretized, respectively, first-order implicit and 
The instantaneous drag and lift coefficients Cp and Ci are second-order central schemes. Equation (8) governing the oscillatory 
defined, respectively, from the drag force Fy and lift force Fy as motion of the square cylinder integrated using fourth-order 
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Runge-Kutta method implemented via a user-defined function (UDF) 
within Ansys Fluent environment. The unsteady solution was sought 
for the fluid subjected to incoming velocity and temperature, and the 
initial cylinder was positioned within the dynamic mesh region. The 
hydrodynamic forces were computed at each time step. The net y- 
direction hydrodynamic force was evaluated to determine the cylinder 
displacement and velocity at that instant time step (Sen and Mittal, 
2011; Willden and Graham, 2006). The cylinder displacement was 
used to adjust the cylinder to the new position, while the velocity was 
accounted as the no-slip cylinder surface velocity in the successive 
time step solution of the fluid equations [Eqs. (5)-(7)]. The time 
advancement was continued by the above process until a quasi-steady/ 
periodic response was achieved. The Courant-Friedrichs-Lewy (CFL) 
number is 0.5 for the entire computational domain (Kumar and Sen, 
2021). For setting the convergence criterion, residual value sets to be 
le-6 for continuity, x-velocity, y-velocity, and energy. 


D. Mesh sensitivity and validation 


Three different meshes with the same O-grid type meshing scheme 
and finer mesh near the cylinder’s surfaces are implemented. Crucial 
parameters such as vibration amplitude A/D, average drag coefficient 
Cp, and average Nusselt number Nuaꝶ are tabulated in Table II. Mesh- 
I, Il, Ill, and IV have, respectively, 10 344, 50 234, 72453, and 80234 
nodes. The difference in the results obtained from Mesh II (baseline), 
and the finer Mesh III is small. The time increment leading to the spatial 
resolution is kept at At = 0.01, and solutions obtained from Mesh II 
are mesh independent and will be used in the subsequent study (with 
slight modifications to accommodate the cylinders’ filleted corners). 

Percentage errors show the deviation of generated values from 
Sun et al. (2021). The present predictions [Eig. 2(b)] at Re= 150, 
m* = 2, and Ur = 10 are compared and in close agreement with 
those of Sun ef al. (2021). Standard deviation from the actual values 
has been plotted in Fig. 2(a). The highest uncertainty can be observed 
in the case of vibrational amplitude, and it depends on the mesh 
refinement, time step, and iteration per time step. Mesh sensitivity 
brings robustness to the numerical simulations and can predict the 
uncertainties in FIV (Wang et al., 2018). 


III. RESULTS AND DISCUSSION 


Both cylinders are elastically supported with m* = 10, and r* is 
changed from 0 (square) to 1 (circular cylinder). For higher amplitude, 
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the structural damping constant is set to be zero (č). Forced convec- 
tion is carried out over the cylinders at a constant Re=150 and 
Pr = 0.7. Reduced velocity (Ur) is varied from 2 to 10 by changing the 
value of spring constant k. Two spacing ratios (L/D) = 2and4 are 
investigated for the detailed analysis of flow separation, recirculation 
behind cylinder, and vortex shedding with three distinctive flow 
regimes. The effect of F on flow characteristics is demonstrated by 
vibration amplitudes, flow contours, streamlines, f/f, frequency 
response, and phase diagrams. 


A. Flow topology for L/D =4 


Figure 3 shows the effect of corner radii with changes on aerody- 
namics and heat transfer parameter. A pattern can be observed for all 
values of A/D in Fig. 3(a); that is, there is the minimal vibrational 
amplitude at lower values at r*=1 and 0.75 for low Ur. At Ur = 6, for 
r* = and 0.75, maxima is achieved for Uc and De (Khan et al., 
2022). Vibrational amplitude of Dc is higher than Uc for both of these 
cases. An interesting phenomenon observed for r* = 0 and 0.5 
between Ur = 6 ~ 8. Maximum values of A/D have shifted from 
Ur = 6 to Ur = 8 for r* = 0 and 0.5 (Zhao and Zhao, 2019). From 
the graph, it is observed that for a particular Ur, the vibrational ampli- 
tude reaches to its maximum value and after it dips to a lower value. A 
generalized assertion can be deduced that rigorous vibration occurs in 
downstream cylinder as compared to upstream cylinder. This phe- 
nomenon is known as wake-induced vibration (Bao ef al., 2012). 

Figure 3(c) shows the Cp variation with Ur for different corner 
radii. Highest Cp is observed for square cylinder and having bigger 
wake at rear of the Uc. For Dc, lower Ur generates negative value of 
average drag coefficient (Cp) as cylinders are not oscillating (Adeeb 
et al., 2018). There is a low-pressure region behind the upstream cylin- 
der, which generates a suction for downstream cylinder. A net force is 
generated from right to left on downstream cylinder, as the low- 
pressure region behind the downstream cylinder has lower vibrational 
amplitude as compared to upstream cylinder. For r* = 0-1 & Ur = 2 
~ 4, there is huge difference in Cp for Uc and De. As Ur increases 
from 4 to 6, lock-in occurs for r = 0.75 and 1 causing Uc to vibrate at 
a certain amplitude. As the vortex sheds from the Uc hits Dc causing it 
to vibrate more rapidly than Dc. A sudden peak in the value of Cp for 
Do can be noticed. Due to the shape change, lock-in for ? = 0 and 0.5 
occurs around Ur = 8, a sharp change can be seen in Cp from Ur = 6 
to Ur = 8 for downstream cylinder. The fluctuation of downstream 


TABLE Il. Numerical simulation results for Re = 150, m* = 2, and Ur = 10. FIV parameters such as maximum amplitude, average drag coefficient, and Nusselt Number are 


summarized. 
Mesh Nodes A/D (Uc) A/D (Dc) Cp(Uc) Cp (De) Nia (Lc) Nuayg(De) 
Mesh-I (coarse) 10344 0.113 0.619 1.102 1.067 5.424 9.359 
Error percentage 12.1% 4.29% 4.29% 5.65% 6.78% 5.93% 
Mesh-II (baseline) 50234 0.125 0.647 1.150 1.096 5.611 9.632 
Error percentage 3.0% 2.8% 5.66% 3.09% 3.57% 3.18% 
Mesh-III (fine-I) 72453 0.127 0.648 1.170 1.125 5.712 9.723 
Error percentage 1.55% 0.14% 4.01% 0.53% 1.8% 2.27% 
Mesh-IV (fine-II) 80234 0.127 0.648 1.170 1.125 5.712 9.723 
Error percentage 1.55% 0.14% 4.01% 0.53% 1.8% 2.27% 
Sun et al. (2021) 15999 0.129 0.6490 1.219 1.131 5.819 9.949 
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FIG. 2. (a) Standard deviation for the computational measurement for four different mesh sizes at various parameters. (b) Comparison of results with publish results of Sun 


et al. (2021) at for Re = 150, m* = 2, and Ur = 10. 


cylinder’s Cp is more due to the function of upstream wake (Alam 
et al., 2003). 

Figure 3(b) shows the effect of r* on Cz with varying reduced 
velocity. Vibrational amplitude and C; are correlated, and the effect of 
it can be seen in given plots. At Ur = 2 ~ 4, there is a very minimal 
coefficient of lift for r* = 0-1 except some of the cases. Downstream 
cylinder of r* = 0 and 1 has higher CI as both cases acting as a single 
body (reattachment) with weak vortex shedding past downstream cyl- 
inder. At Ur = 6, for r* = 0.75 and 1, vibrational amplitude A/D and 
lift coefficient Cz shoot up as lock-in occurs. Cz is higher for the 
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downstream cylinder because vortex sheds from the Uc impinge on 
the front surface of the Dc. T = 0.75 breaks the steady flow regime to 
enter in unsteady flow. 

Figure 4 shows different Z-vorticity contours for all the defined 
cases at different Ur. Detachment of shear layer followed by rolling of 
those shear layer due to wake behind the cylinder causes vortex shed- 
ding. For square cylinders at Ur = 2 ~ 6, the shear layer separated 
from Uc gets reattached to Dc causing it to behave like a single cylin- 
der with unsteady flow (Yao and Jaiman, 2019). Separation occurs 
from leading edge of Uc attached again to leading edge of Dc. 
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FIG. 3. Effect of f on FIV of the tandem cylinder for L/D = 4 by varying Ur (a) A/D vs Ur, (b) Ci vs Ur, (c) Cp vs Ur, and (d) Nuavg vs Ur. 


Two square cylinders are not vortices in alternate manner with inline 
formation. This exhibits the 2S vortex shedding from two cylinders 
together. For Ur = 8 and r* = 0, shear layers are pushed into two par- 
allel streets by the vibration of downstream cylinder causing C(2S) 
past De. 

Figure 5 shows the streamlined structure of wakes behind the cyl- 
inders. Streamlines for square cylinder show a varying wake for all 
reduced velocities behind the Uc. For a square cylinder, at higher Ur, 
wake behind the Uc becomes significantly low as vortices from Uc 
start dominating the region. Even in reattachment regime, there are 
some distinctions. As it can be observed for r = 0.5 and 0.75 at 
Ur = 2 ~ 6, flow getting separated at the trailing edge of the Uc and 
getting reattached at Dc. Steady flow can be observed past down- 
stream cylinder. As shown in Fig. 5, wake at r* = 0.5 at Ur = 6 is sym- 
metrical for both Uc and De. 

Flow characteristics of this set are similar to the steady flow past 
an isolated cylinder. No vortex shedding is observed past downstream 
cylinder. Flow topology becomes different at r* = 0.75, and higher 
vibrational amplitude of Do becomes massively high. In the case of cir- 
cular cylinder, tandem cylinders at Ur = 2 and 4 act as single cylinders 
due to reattachment, and 2S vortex shedding can be observed for the 
Ur = 2. r = 1 and Ur = 4 vortices from the upper surface of cylinder 
and lower surface of Dc. A major similarity is observed for ? = 0 and 
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1 at Ur = 2, and wake structure is changing behind upstream cylinder, 
but no actual presence of wake is observed at the downstream cylinder. 

For r = 1 and Ur = 6, physics of flow becomes more dynamic 
as the cylinder is their peak of vibrational amplitude as shown in 
Fig. 3(a). Vortices from Uc are impinging on the Dc periodically, and 
vortices from upper surface of Uc interact with separation at the lower 
surface of Do. This interaction causes vortex stretching nearby the 
downstream cylinder and lies under “co-shedding” (Griffith et al., 
2017). Vortex shedding from tandem cylinders changes to 2S or Von 
Karman vortex shedding. As Ur reaches the threshold of 6, unsteady 
flow starting to emerge for r = 0.75 and 1 similarly results started 
showing up for r = 0 and 0.5 at Ur = 8. For square cylinder past 
Ur = 8, unsteady flow is actively present. As show in Fig. 4, vortex 
shedding from cylinders start as 2S, and in later stages, it turns to be 
C(2S). There is reasonable gap in vortices as De is currently under 
lock-in period. For r* = 0.5, it does occur at Ur = 8, as compared to 
square cylinder vortex shedding amplitude, which is significantly low. 
This causes vortices to shed in an inline manner, and vortices are 
formed in combination of C(2S) and 2S. A sharp drop is observed for 
all cases after the lock-in period. Most of the cases from r* = 0.5 ~ at 
Ur = 8 ~ 10 shows similar flow structure behind the De. 

Figure 6 shows the instantaneous velocity and pressure contours 
for r = 0-1. In the case of square cylinder at Ur = 6, there is a 
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FIG. 4. Z-vorticity contours for varying ? = 0-1 and Ur = 2 — 10 for L/D = 4. 


prominent low pressure between the cylinder, which can be confirmed 
by negative value of Cp. As fluctuation starts in upstream cylinder 
causes a varying low-pressure region between the cylinders and past 
downstream cylinder another fluctuating wake can be observed which 
is similar to Yao and Jaiman (2019). Although when Uc and Dc both 
vibrate, it causes a larger wake as compared to other cases. This differ- 
ence is prominent in Ur = 6 and Ur = 8 for r* = 0. Other configura- 
tions are similar in nature, but values fluctuate on the bases of type of 
vibration as in the case of Ur = 6 for r = 0.5 and 0.75, and velocity 
and streamlines show that steady flow does not fluctuate more as com- 
pared to crossheading. Although for circular cylinder at Ur = 6, co- 
shedding phenomenon is taking place; hence, a huge wake is formed 
behind the downstream cylinder. 

Vortex shedding imprints are captured by contours of fs /fn and 
power spectral density of lift coefficient for L/D = 4. Figure 7 shows 
the ratio of natural frequency and vortex shedding frequency as dotted 
line shows the lock-in period where ff = 1. From the Figs. 3(a) and 
7(a), it is evident that lock-in for Uc and De occurs at Ur = 6 
(Williamson, 1996). Synchronization remained persisted at the same 
values of Ur for r* = 0.75. As corner radii stared to decreases, there is 
shift start occurring for r* = 0.75 and 1 as dotted line tilts toward the 
Ur = 8. Downstream cylinders try to resist the “shift,” but fluctuation 
in Uc will cause changes in Dc. 
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Figures 8(a)-8(d) illustrate the shedding frequency for r* = 0-1 
at Ur = 6. Shedding frequency of r* = 0 is less than the natural fre- 
quency (fa = 0.166). r* = 0.5 generates identical results as r* = 0. 
Values of downstream cylinder and upstream cylinder frequency are 
almost identical. Shedding frequency at r* = 1 and 0.75 escalates 
quickly to higher values. f* for r* = and 0.75 approaching the natu- 
ral frequency (% = 0.166) A second peak is observed for r = and 
r* = 0.75 illustrating frequency modulation (Fed). Figures 8(e)-8(h) 
depict the shedding frequency at Ur = 8 for r* = 0-1. T = 0 (square 
cylinder) and r* = 0.75 registered the growth in f* while changing 
the reduced velocity. Primary vortex shedding approaches the natural 
frequency (fa = 0.125) with the modulation. At r* = 0, change in 
vortex shedding for Dc is higher as compared to Uc (Borazjani and 
Sotiropoulos, 2009; Prasanth and Mittal, 2009). Vortex shedding fre- 
quency of r* = 0.75 and r* = 1 (Dc/Uc) has surpassed the natural 
frequency (fa = 0.125). Frequency modulation has been omitted 
from circular cylinder at Ur = 8. From Fig. 8, conclusion can be 
drawn on the lock-in phenomenon as lock-in for r* = 0 and 0.5 
occurs as Ur = 8, while for r* = 0.75 and 1 occurs as Ur = 6, 
respectively. 

Figure 9 shows phase diagrams for different corner radii at three 
Ur for Uc and Dc. At Ur = 4, for all cases Cr, and Y* appear to be in 
same phase. At Ur = 6, lock-in occurs, which is also shown in the 
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FIG. 5. Velocity streamlines for different cases showcasing variants of wake behind the cylinders. 


graphs of Uc and Dc at r* = 0.75 and 1. For r* = 1 and 0.75, phase 
lag can be observed at Ur = 8. 

Amplitude for upstream and downstream cylinder has dropped to 
a certain lower value contrary to the nature of Cz, which is stringent to 
higher values. For all cases with Uc and Dc, phase change is occurred 
when Ur changes from 6 to 8 (Lin et al., 2014). Equation (13) will pro- 
vide the direction of force, which dictates the position of cylinder 


Fy [ANV 
wo! p (13) 


Figures 10(a)-10(d) show the dependence of Nu on @ for Ur = 6 
at r = 0-1. For square cylinder, the frontal area has higher Nu 
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because frontal area is in direct contact of upcoming flow. At 0 = 45°, 
it depicts rise for Nu due to sharp edge. Heat transfer after 45° started 
dropping due to side bubble for r* = 0. Low Nu at the frontal part of 
Dg (r* = 0) due to the presence of “dead water” between the cylinders. 
At r = 0.75, a smooth curved is followed as started from lower to 
higher Nu at corner of geometry. Nu dropped to lower values after 
0 = 135° due to shear layer separation, which reattached to down- 
stream cylinder around 0 = 35° — 45°. As shown in Fig. 10, separa- 
tion from downstream cylinder is almost at the same angle as 
upstream as separation is dependent on the geometry. 

Almost same characteristics are shown by r* = 0.5 and r* = 0.75, 
as only flow separation angle has been shifted from 135° to 120°. 


35, 013607-10 


Physics of Fluids ARTICLE scitation.org/journal/phf 


r*=0 r*=0.5 r*=0.75 r*=] 
| | | 


rye 20 2 ͤ—ümä 
M w a 
Ur=6 ae * ae og tf og 
M P A 
5 > a 5 71>, ar 5 ö 
S 


S 
Pressure 
Co 


Ur=10 


wo 


— 


* 


Velocity Magnitude: 0.00 0.17 0.33 0.50 0.67 0.83 1.00 1.17 1.33 1.50 
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For circular cylinder, a perfect curve is followed as 0° have highest Nu. Figures 10(e)-10(h) depict the polar graph between Nu and 6 for Ur 
Along the circular cylinder, sudden drop can be observed in Nu for cir- = 8 at r = 1. The purpose of illustrating polar graphs at two differ- 
cular cylinder indicating flow separation due to adverse pressure gradi- ent reduced velocities is to understand the effect of vibrational 
ent. Separation angle for Uc is around 125°, while for De, it is at 135°. amplitude at transitional level. For r = 0, Nu(Dc) have increased 
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FIG. 9. Phase diagram of upstream cylinder (blue grid) and downstream cylinder (red grid) at Ur = 2, 4, and 8. 


due to vibrational amplitude although position of point of separation movement synchronization with upstream cylinder. For circular cyl- 
has not changed. Core of r* = 0.5 and r* = 0.75 for downstream cyl- inder at Ur = 8, vibrations have suppressed for upstream cylinder in 
inder has changed as downstream, which has shown vibration. At contrast to downstream, which leads to the interaction of shear 
Ur = 6 and r* = 0.75 have sufficient but downstream cylinder is in layers from Uc with Dc. 
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Figure 11 shows the non-dimensional temperature variation 
around the cylinder when it is maintained at a constant temperature. 
Upcoming flow at T* = 0, Re = 150, and Pr = 0.7 and cylinders’ sur- 
face are maintained at T* = 1. White color (T* = 0) in domain signi- 
fies the presence of unaffected area from heat transfer around the 
cylinder. Nuyg of downstream cylinder is highly dependent on the 
Nuayg of upstream cylinder (Alam et al., 2021; 2020). Overall, Nuayg is 
dependent on flow topology and vibrational amplitude. Higher vibra- 
tional amplitude led to the higher Nuag. Figure 3(d) shows the 
Nusselt number dependence on Ur and H. In the case of square cylin- 
der, Nuayg is relatedly low as compared to other cases; at ? = 0, flow 
separation occurs at the trailing edge of square cylinder, which is early 
separation that causes poor heat transfer (Dhiman et al, 2006). 
Heated flow from the upstream cylinder reattaches to downstream cyl- 
inder, causing poor heat transfer from De for Ur = 2 ~ 4 (Zhou and 
Yiu, 2006). Even for Do, delayed flow separation causes better heat 
transfer (Dwivedi and Dhiman, 2019). 

As soon as cylinder reaches in lock-in period, a rise can be seen 
in Nuaw for Uc and Dg. The rise in Dc is more evident and sharper 
because Do is getting exposed to fresh flow, which led to better heat 
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transfer. At r = 0.75~1 for Ur = 6 and r* = 0.75~1 for Ur = 8 
and, a sharp change in Nua can be observed for Dc. Same shift 
can be observed in the value of Nua. Vibrational amplitude of Dc at 
r* = at Ur = 8 is significantly higher than r = 0.50 at Ur = 8, but 
Nuayg of r* = 0 is slightly lower than r* = 0.50 at Ur = 8. This defines 
the flow topology over the downstream cylinder and deifies the com- 
pleted dependence of Nia at A/D (Sun et al, 2021). 


B. Flow topology for L/D = 2 


Figure 12 shows corner radii effect on the FIV and heat transfer 
parameters at L/D = 2. A general trend can be observed for the 
upstream and downstream cylinder in Fig. 12(a); that is, a very low 
vibrational amplitude is observed at Ur = 2 ~ 4 for * = 0~1. As 
enormous rise in A/D can be seen for Uc and Dc at r = 1, vortex 
shedding frequency approaches the vortex shedding frequency 
(Khalak and Williamson, 1997). 

At Ur = 6, maximum value of Uc and Dg is achieved for circular 
cylinder. A significant change has been observed in r* = 0.5 and 0.75 
as vortices from Uc start hitting Dc, although due to less spacing 
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FIG. 11. Non-dimensional temperature contour for both heated cylinders at a constant wall temperature with varying Ur and F. 
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FIG. 12. Effect of r* on FIV of tandem cylinder by varying Ur. (a) 4 vs Ur, (b) Ci vs Ur, (c) Cp vs Ur, and (d) Mag vs Ur at L/D = 2. 


vortices are not able to roll up properly (Lin et al., 2014). A/D of De is 
usually higher than the A/D of Uc. For square cylinder, no vibration 
is observed for Ur = 2 ~ 6. At Ur = 8, r* = 0.5 and 0.75 achieve its 
maximum value, which is same as r* = 1 at Ur = 6. A “shift” is 
observed in values as corner radii are changed from 1 to 0 (Zhao and 
Zhao, 2019). For square cylinder, maxima is observed at Ur = 8, but 
value is comparatively low as compared to other geometries. At 
Ur = 10, A/D is plummeted to lower values for all cases. This con- 
firms usual behavior of A/D with respect to Ur. A shift is caused for 
peaks of A/D, which is produced due to change in corner radii. For 
some cases, a nominal difference in vibrational amplitude can be 
noticed between Uc and De. Vortex shedding from upstream cylinder 
impinges on the downstream cylinder, which causes fluctuation in all 
parameters (Borazjani and Sotiropoulos, 2009). 

Figure 12(c) illustrates the Cp with respect to Ur at varying cor- 
ner radii. Cp, in general, is function of forces applied on the body in 
axial direction. Cp of Dc is also dependent on the Uc. Maximum Cp 
is observed for upstream square cylinder at Ur = 10, while minimum 
Cp is observed for Ur = 2 and 4. Negative values of Cp can observed 
for downstream cylinders at lower Urs, as Dc captured in wake of Uc. 
Existence of “dead water” between the cylinders shows very low- 
pressure region, which is very low in comparison with wake of Dc. 
Net force opposite in axial direction is experienced by the De. After 
Ur = 6, Dc is subjected to vibration consequently Cp started rising. 
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Figure 12(b) shows the changes in C; with respect to reduced velocity 
for different cases. Vibrational amplitude and Cz are concurrent, and 
the effect of it is validated from above graphs. At Ur = 6, peak in Ci is 
observed for all cases. Maximum of C is observed at Ur = 6 for De, 
while lowest is observed for Ur = 2 for r* = 1. 

Figure 13 shows the time history of A/D for r = 0-1 in lock-in 
phase for downstream and upstream cylinder. In most of all the cases, 
amplitude of downstream cylinder is higher as compared to upstream 
cylinder. Fluctuations in shear layer of upstream cylinder also affect 
the vibrational modes of downstream cylinder. A/D for r = 1 is 
peaked at Ur = 6 as reduced velocity is increased to Ur = 8, and 
amplitude drops to lower values, but CA persists its nature. 

Figure 14 depicts different Z-vorticity contours for all the defined 
cases at different Ur. Detachment of shear layer followed by rolling of 
those shear layer due to wake behind the cylinder causes vortex shed- 
ding. For r* = 0, Uc overshoots the vortices over Dc and hence acts as 
“single body,” which is similar to isolated cylinder. 2S type of vortex 
shedding is observed in these particular cases. Although this regime is 
dependent on the vibrational modes of both cylinders, r = 0.5, 0.75, 
and 1 also show single body regime at lower values of vortex shedding. 
For r = 0.5 at Ur = 6, as amplitude of Dc overshoots due to the 
lock-in, which changes the vortices structure in the wake. Shear layer 
separated from upstream layer is pushed in opposite directions (Sun 
et al., 2021). In the case of circular cylinder, Dc pushes the vortices to 
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Ur = 6. 


the extreme and causes a parallel vortex street. This vortex street is 
also known as C(2S) vortex shedding (Williamson and Roshko, 1988); 
hence, De also dictates the nature of vortex shedding. As A/D of Dc 
plummeted with increase in Urs, vortex shedding is similar to isolated 
cylinder. A phenomenon of vortex stretching can be observed for r* 
= 0.5 and 1. 

Figure 15 shows the velocity streamlines for some selected cases 
with varying Ur. For square cylinder at low, Ur shows symmetric 
wake between the cylinders, while at Ur = 10 does not illustrate a sta- 
ble wake due to the vibration of De. A side bubble can be noticed for 
square cylinder, which caused due to edge separation. r* = 0.5 shows 
the symmetric wake (recirculation) between cylinder although due to 
the Coanda effect, flow follows the topology. Topology of figure delays 
the separation, and side bubble is not visible in this configuration. 
Oscillating circulation is observed at the rear of Dc. For r* = 0.75 and 
Ur = 4, it shows asymmetric circulation between upstream and down- 
stream cylinder. Higher vibration amplitude destabilizes the wake and 
eliminates the low-pressure region between the cylinders. Circular cyl- 
inder has highest of vibrations for downstream cylinder, and a promi- 
nent difference can be observed at Ur = 6 and Ur = 8. Figure 15 
illustrates that vortex for Ur = 6 is pushed in a parallel manner, while 
vortices for Ur = 8 are in inline manner. 

Figure 16 instantaneous pressure and velocity contour plots for 
r = 0-1 at L/D = 2. Spacing ratio is reduced to 2; hence, mostly sin- 
gle body effect can be noticed in most of the cases. “Dead water” can 
be observed between the cylinders when Dc is not having dynamic 
response. Due to 2S type of vortex shedding, the pressure behind the 
downstream cylinder is fluctuating. r* = 1 and Ur = 6 have larger 
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wakes due to C(2S) type of vortex shedding. Downstream cylinder 
pushes the shear layers from upstream cylinder to create two parallel 
streets; hence, low-pressure region is created between the parallel 
streets. 

Vortex shedding characteristics of vortex shedding are captured 
by contours of f/f, and power spectral density in Figs. 17 and 18. 
Figures 17(a) and 17(b) show the ratio of natural frequency and vortex 
shedding frequency. Lock-in period exists around f;/f;, = 1. From the 
Figs. 17(a) and 17(b), it is evident that lock-in for Uc and De occurs 
around Ur = 6and8 (Mysa et al, 2016). Synchronization remained 
persisted at the same values of Ur for r* = 1, 0.50, and 0.75. As corner 
radius starts to decrease, the shift does not occur in case for L/D = 2 
as compared to/D = 4. Both upstream and downstream cylinders try 
to resist this shift, but fluctuation in Uc will cause change in Dc. 

Figures 18(a)~18(d) illustrate the primary shedding frequency for 
r* = 0-1 at Ur = 6(L/D = 2). Shedding frequency of r = 0 is lower 
than the natural frequency (fa = 0.166) although values have 
increased as compared to L/D = 4. r* = 0.5 results differ from as r* 
= 0, f*(r* = 0.5) escalated quickly and approaching the n showing a 
huge leap in vibration amplitude as well. Values of downstream cylin- 
der and upstream cylinder frequency are similar in terms of absolute 
values. 

Shedding frequency at r* = 1 and 0.75 shows similar values at 
Ur=6. f* for r* = 1, 0.75, and 0.5 impending the natural frequency 
(fn = 0.166). A second peak is observed for r = 1 and r? = 0.5 illus- 
trating frequency modulation (fmoa) around 0.4. Figures 18(e)-18(h) 
depict the shedding frequency at Ur = 8 for r* = 0-1. f* =0.125(Uc/ 
De) of r = 0 (square cylinder) exactly overlaps the natural frequency 
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FIG. 14. Z-vorticity contours for varying ? = 0-1 and Ur = 2 — 10 for L/D = 2. 


(fa = 0.125) resulting huge vibrational amplitude. No secondary vortex 
shedding is observed for r = 0 at Ur = 8 (Willden and Graham, 2006; 
Williamson, 1996). Vortex shedding frequency of ? = 0.5, 0.75, and 1 
Oe / Uc) has dropped as compared to the r = 0, which is contrary to 
L/D =4. Frequency modulation can be noticed for r = 0.75 and 1 as 
two extra peaks have been observed. Lock-in for r* = 0 occurs as 
Ur = 8, while for r = 0.75 and 1, it occurs as Ur = 6. At r = 0.5, it 
shows larger range for lock-in from Ur = 6 to Ur = 8. 

Figure 19(a) shows phase diagrams for different corner radii at 
three Ur for Uc and De. At Ur = 4, for all cases Ci and Y* appear to 
be in same phase. At Ur = 6, lock-in occurs, which is also shown in 
the graphs of Uc and De at r* = 0.75 and 1. Amplitude for upstream 
and downstream cylinder has dropped to a certain lower value con- 
trary to the nature, which is stringent to higher values (Papaioannou 
et al., 2008). Complete frequency synchronization occurs at Ur = 8 in 
which the vortex shedding frequency collapses exactly upon the natu- 
ral frequency of the vacuum; hence, the cylinder moves along a peri- 
odic orbit in Fig. 19(a). This occurred for Uc and Dc at Ur = 8, and 
“eight” type figures can be observed for r* = 0.5,0.75, and 1. For r* 
= 0.5, Uc phase change occurs at Ur = 6 contrary to other configura- 
tions of Uc and De, which have opposite direction of amplitude and 
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lift due to phase change. Figure 19(b) shows the instantaneous position 
of cylinder. It depicts the positions of Uc and Do, which are opposite 
sign when one cylinder is at maxima. 

The Nusselt number is dependent on the spacing ratio, corner 
radii, and flow topology around the cylinders (Zhou and Yiu, 2006). 
Figures 20(a)-20(d) demonstrate the dependence of Nu on 0 for 
Ur = 6 at r = 0-1 (L/D = 2). Spacing ratio does not affect the Nu 
for the frontal area of cylinder as square cylinder has higher Nu, which 
is same as L/D = 4. At 0 = 45°, it depicts the rise for Nu due to sharp 
edge. Heat transfer after 45° started plummeting due to side bubble 
for r = 0. Rear separation can also be observed for square cylinder 
around 0 = 135°. For r = 0.5, Dc has registered a growth in Nu as 
compared to the L/D = 4 due to shift in vibrational amplitude. 
Nu(D¢) for r* = 0.5 started dropping after 0 = 65° due to shear layer 
separation. As shown in Fig. 20, separation from downstream cylinder 
is almost at the same angle as upstream as separation is dependent on 
the geometry. For spacing ratio of 2, angle of reattachment and angle 
of separation are entirely different due to overshooting of vortices 
from Uc. Hence, spacing ratio dictates the separation point from 
downstream cylinder. For r* = 0.5, separation at L/D = 4 occurs 
around 0 = 45° — 50°, while for L/D=2, it has shifted to 
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FIG. 15. Velocity streamlines for different cases show casing variants of wake behind the cylinders for spacing ratio of 2. 


@ = 65° — 70°. Identical characteristics are shown by r* = 0.5 and r* 
= 0.75. For circular cylinder, a perfect curve is followed as 0°, which 
has highest Nu. Along the circular cylinder, sudden drop can be 
observed in Nu for circular cylinder indicating flow separation due to 
adverse pressure gradient. 

Figures 20(e)-20(h) depict the polar graph between Nu and for 
Ur = 8 at r* = 0-1. The purpose of illustrating polar graphs at two 
different reduced velocities with spacing ratio is to comprehend the 
effect of vibrational amplitude at transitional level. For r* = 0, Nu(Dc) 
has three points of separation for downstream cylinder starting from 
leading-edge separation to trailing edge separation. Point of separation 
has delayed for r* = 0.5 (Dc) and r* = 0.75 (Dc) as compared to 
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upstream cylinder for same configurations. Due to less interaction at 
frontal surface of De (r* = 0.5), it causes drop in heat transfer. At fil- 
leted area (0 = 45°), a rise of Nu is observed generating “butterfly” 
shape in Fig. 20(f). Figure 20 shows the temperature variation around 
the cylinder when it is maintained on a constant temperature. Na of 
downstream cylinder is highly dependent on the Nua of upstream 
cylinder. Heated flow from upstream cylinder interacts with down- 
stream cylinder. Overall, Nua is dependent on flow topology and 
vibrational amplitude (Zhou and Mahbub Alam, 2016). Higher vibra- 
tional amplitude led to the higher Nuayg. Figure 1 1(d) shows variation 
in the Nusselt number with respect to Ur with corner radii. Upstream 
cylinder of square cylinder has least Uc in comparison with other 
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FIG. 19. (a) Phase diagram of upstream cylinder (Blue Grid) and downstream cylinder (Red Grid) at Ur = 2, 4, and 8 for L/D = 2 and (b) time history for amplitude of 


r’=0.5 and Ur = 6. 


configurations. Flow separation for T = 0 occurs at the leading edge 
of square cylinder, which is main cause of less Nilag. Figure 21 dem- 
onstrates the non-dimensional temperature changes around the 
upstream cylinder when it is maintained at a constant temperature. 
Upcoming flow sets to be at T* = 0, Re = 150, and Pr=0.7, while 
both cylinders’ surface are held at T* = 1. White color (T* = 0) in 
domain signifies the presence of unaffected area from heat transfer 
around the cylinder, while affected area is shown by colored contour 
given in Fig. 21. Temperature contours are similar to vortex shedding 
patterns as vortices carry heat after shear layer separation. 

Nuag of Dc is less than Nugyy of Uc due to plummeted flow 
interaction at the rear of upstream cylinder. Circular cylinder shows 
maximum values of Nua due to delayed flow separation causing 
higher heat transfer to shear layers (Zhang et al., 2019). At Ur = 6, r* 
= 1 shows higher values as compared to Ur = 4 for Dc; A/D of Dc 
gained momentum and start interacting to upcoming flow rather than 
flow passing from Uc. Nua of r* = 1 experience maxima at Ur = 8, 
while a sharp drop is observed at Ur = 10. At L/D=2, Nuaꝶ have a 
lower value of Dc as compared to Dc of L/D=4. Better heat transfer 
occurs in “Reattachment” as compared to “single body.” This is defin- 
ing the flow topology over the downstream cylinder and deifies the 
completed dependence of Nua at A/D. 
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IV. CONCLUSION 


The effect of spacing ratio and corner radius on forced convec- 
tion from tandem cylinders is numerically investigated. For all cases, 
major parameters such as Reynolds number = 150, Prandtl number 
=0.7, mass ratio = 10, and damping ratio = O are fixed. Reduced 
velocity varied Ur=2-10 while two different spacing ratio L/D 
= 2and4. Four corner radii (r*) of cylinder varying from 0, 0.5, and 
0.75 to 1 are employed. Major observations are summarized as follows: 


Effect of corner radius illustrates the varying flow topology over 
the cylinders. Square cylinders show the leading-edge separation 
with side bubble, while circular cylinder depicts the separation 
past front surface with primary wake. Delayed flow separation 
can be observed in r* = 1 and 0.75. Symmetric and asymmetric 
wakes have been observed with change in Ur, T, and Z, 

For L/D = 4, maxima of vibrational amplitude are observed for 
De at r* = 1 while minima occurred at Ur = 2 for r* = 0. 
Maximum Cp is noted for r* = 1 for upstream cylinder, which 
occurred due to large wake behind the cylinder. Negative value of 
Cp is discovered, which explains the presence of “dead water” 
between the cylinders. Lock-in for r* = 0 and 0.5 occurs at 
around Ur = 8, while for r* = 0.75 and 1, it occurs at Ur = 6. 
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FIG. 21. Temperature contour for both heated cylinder at a constant wall temperature with varying Ur and F. 


This phenomenon is termed as “shift.” Direction of forces on cyl- 
inder changes due to change in Ur. At Ur = 4, for all cases Ci. 
and Y* appear to be in same phase. At Ur = 6, lock-in occurs, 
which is also shown in the phase diagrams of Uc and De. At r* 
= 0.75 and 1 (Ur=6) are in same phase while phase lag 
observed at Ur = 8. For all cases, phase change is occurred when 
Ur changes from 6 to 8. 

Different type of vortex shedding has been observed along with 
various flow patterns. Flow patterns are closely related to shear 
layer separation. Changes in shear layer separation can cause com- 
plicated wake interaction and chaotic shedding. Except square 
cases, most of the cases have shown reattachment after separation. 
Vortex shedding from tandem cylinder is either C(2S) or 28. 

For Ur = 6, F for r = 1 and 0.75 approaching the natural fre- 
quency (% = 0.166). A second peak is observed for r = 1 and r* 
= 0.75 illustrating frequency modulation. At Ur = 8, vortex shed- 
ding frequency of r = 0.75 and r* = 1 (Dc/Uc) is impending the 
natural frequency (fa = 0.125) causing synchronization. Frequency 
modulation can be observed from the power spectral density. 

For square cylinder, the frontal area has higher Nu because fron- 
tal area is in direct contact of upcoming flow. At r = 0.75, a 
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smooth curved is followed as started from lower to higher Nu at 
corner of geometry. For circular cylinder, a perfect curve is followed 
as 0° have highest Nu. Along the circular cylinder, sudden drop 
can be observed in Nu for circular cylinder indicating flow separa- 
tion due to adverse pressure gradient. Vibration amplitude affects 
the local Nu due to change in vibrational modes. Forced convection 
on the tandem cylinders in this study is dependent on major factors 
such as A/D, flow topology, spacing ratio, and corner radii. Higher 
vibration in tandem cylinders adds in the value Nua. Maximum 
value of Nua is at r = 1 and Ur = 6 for L/D = 4 (Uc), while 
minimum is recorded at r* = 0 and Ur = 2 (De). 

At L/D = 2, flow characteristics becomes simpler. Maximum A/D 
is observed for r* = 1 (Dc) at Ur = 6. For r = 1 and Ur = 6, lock- 
in occurs along with C(2S) type of vortex shedding. Phases illustrate 
the force’s direction and force response on vibrational amplitude. At 
Ur = 4, Ci and Y* appear to be in same phase. Complete frequency 
synchronization occurs at Ur = 8 in which the vortex shedding fre- 
quency collapses exactly upon the natural frequency of the cylinder; 
hence, the cylinder moves along a periodic orbit. 

Lock in period exist around f,/f,, = 1. Contour graphs supported 
that phenomenon (lock-in) for Uc and De occurs around 
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Ur = Gand 8. Synchronization remained persisted at the same 
values of Ur for r = 1, 0.50, and 0.75. The “shift” did not occur 
in case for L/D = 2 as compared to/D = 4. Both upstream and 
downstream cylinders try to resist this shift, but fluctuation in 
Ug will cause changes in De. 

Steady flow along with “single body” regime can be observed at 
lower Urs. For Ur = 6 at r* = 0.5 and 1, it depicts C(2S) type of 
vortex shedding, which occurred due to interaction of vortices 
from lower surface and upper surface of Dc. A parallel vortex 
street is created by pushing the shear layer, which is approaching 
from upstream cylinder. 

For L/D = 4, separation from downstream cylinder is almost at 
the same angle as upstream as separation is dependent on the 
geometry. For spacing ratio of 2, angle of reattachment and angle 
of separation are entirely different due to overshooting of vortices 
from Uc. Hence, spacing ratio dictates the separation point from 
downstream cylinder. For r = 0.5, separation at L/D = 4 occurs 
around 0 = 45° — 50°, while for L/D = 2, it has shifted to 
0 = 65° — 70°. Along the circular cylinder, sudden drop can be 
observed in Nu for circular cylinder indicating flow separation 
due to adverse pressure gradient. The purpose of illustrating 
polar graphs at two different reduced velocities with spacing ratio 
is to comprehend the effect of vibrational amplitude at transi- 
tional level. For r* = 0, Nu(Dc) have three points of separation 
for downstream cylinder starting from leading-edge separation to 
trailing-edge separation. For L/D = 2, maximum Nu occurs at 
= 1 and Ur = 6 (Uc), while minimum lies on r = 0 and Ur 
= 2 (Dc). At L/D = 2, Nuayg have a lower value of Dc as com- 
pared to Dc of L/D=4. Better heat transfer occurs in 
“reattachment” as compared to “single body.” 
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